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Q_i. Abstract 



Hadronic charmless B decays to scalar mesons are studied within the framework of QCD fac- 
torization (QCDF). Considering two different scenarios for scalar mesons above 1 GeV, we find 
that the data favor the scenario in which the scalars ao(1450) and -K"q(1430) are the lowest lying 
qq bound states. This in turn implies a preferred four-quark nature for light scalars below 1 GeV. 
Assuming ^(1430) being a lowest lying qs state, we show that the data of B — > Kq (1430)t/') and 
B — > Kq (1430)(p, cj, </>) can be accommodated in QCDF without introducing power corrections in- 
duced from penguin annihilation, while the predicted B~ — > Kq (1430)7r~ and B — > Kq (1430)7r + 



are too small compared to experiment. In principle, the data of Kq (1430)7r modes can be explained 
if penguin-annihilation induced power corrections are taken into account. However, this will destroy 
the agreement between theory and experiment for B — > Kq (1430)(?/ \p,uj,(p). Contrary to the 
pseudoscalar meson sector where B — > Krf has the largest rate in 2-body decays of the B meson, 



X 



we show that B{B -)■ KqI]') < B(B -> KqT]). The decay B° -)■ 00(980)+^" is found to have a rate 
much smaller than that of B° — > ao(980) + 7r~ in QCDF, while it is the other way around in pQCD. 
Experimental measurements of these two modes will help discriminate between these two different 
approaches. Assuming 2-quark bound states for /o(980) and /o(500), the observed large rates of 
/o(980)K and f (980)K* modes can be explained in QCDF with the / (980) -/ (500) mixing angle 
9 in the vicinity of 20°. However, this does not necessarily imply that a 4-quark assignment for 
/o(980) is ruled out because of extra diagrams contributing to B — > /o(980)-RT(*) . Irrespective of 
the mixing angle 9, the predicted branching fraction of B° — > fo(980)p° is far below the Belle 
measurement and this needs to be clarified in the future. 
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I. INTRODUCTION 



In the past few years there are some progresses in the study of charmless hadronic B decays 
with scalar mesons in the final state both experimentally and theoretically. On the experimental 
side, measurements of B decays to the scalar mesons such as /q(980), /o(1370), /o(1500), /o(1710), 
ao (980), ao (1450) and ^(1430) have been reported by BaBar and Belle; see Tables HI and [TT1 for 
a summary of the experimental results. It is well known that the identification of scalar mesons 
is difficult experimentally and the underlying structure of scalar mesons is not well established 
theoretically. The experimental measurements of B — > SP and B — > SV, where S, P, V stand 
for scalar, vector and pseudoscalar mesons, respectively, will provide valuable information on the 
nature of the even-parity mesons. On the theoretical side, hadronic B decays to scalar mesons 
lave been studied in the QCD-inspired approaches: QCD factorization (QCDF) and pQCD 
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2C]. 



In this work, we would like to revisit the study of the 2-body charmless decays B — > SP and 
B — > SV within the framework of QCDF for the following reasons: (i) In [6J we have missed some 
factorizable terms (more precisely, the /o and a{] emission terms) in the expressions for the decay 
amplitudes of B — > foK, o,qt:, a^K. (ii) Attention has not been paid to the relative sign difference 
of the vector decay constants between clq and clq and between Kq and Kq or Kq~ and Kq + in our 
previous study, (iii) There were some errors in our previous computer code which may significantly 
affect some of the calculations done before, (iv) Progress has been made in the past in the study of 
B — > S transition form factors in various approaches 0, [21-26]. (v) Experimental data for some 
of B — > SV decays such as Kq (1430)0, Kq (1430)/? and Kq(1A30)uj are now available, (vi) It is 
known that in order to account for the penguin-dominated B — > PP, VP, VV decay modes within 
the framework of QCDF, it is necessary to include power corrections due to penguin annihilation 



271 . [28J. In the present work, we wish to examine if the same effect holds in the scalar meson 



sector; that is, if the penguin-annihilation induced power corrections are also needed to explain the 
penguin dominated B — > SP and B — > SV decays. 

This paper is organized as follows. We specify in Sec. 2 various input parameters for scalar 
mesons, such as decay constants, form factors and light-cone distribution amplitudes. The relevant 
decay amplitudes are briefly discussed in Sec. 3. Results and detailed discussions are presented 
in Sec. 4. Conclusions are given in Sec. 5. We lay out the explicit decay amplitudes of B° — > 
(/o, a[})(ivr, 7r) in Appendix A. 



II. PHYSICAL PROPERTIES OF SCALAR MESONS 

In order to study the hadronic charmless B decays containing a scalar meson in the final state, 
it is necessary to specify the quark content of the scalar meson. For scalar mesons above 1 GeV 
we have explored in 0] two possible scenarios in the QCD sum rule method, depending on whether 
the light scalars Kg (800), ao(980) and /o(980) are treated as the lowest lying qq states or four- 
quark particles: (i) In scenario 1, we treat Kq (800), ao(980), /o(980) as the lowest lying qq states, 
and A"q (1430), ao(1450), /o(1500) as the corresponding first excited states, respectively, and (ii) 
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TABLE I: Experimental branching fractions (in units of 10 ) of B decays to scalar mesons [lj 



Mode 




BaBar 


Belle 


Average 


B(B+ 


->■ / (980)i<:+)B(/o(980) -> tt+tt - ) 


10.3 ±0.5+? ° 

— 1.4 


8.8±0.8t°-^ 


9.4+^ 

— l.U 


B(B° 


/o(980)K°)B(/o(980) tt+tt") 


6.9 ± 0.8 ± 0.6 


7.6 ±1.7+?'? 

' — 1.3 


7.0 ±0.9 


B(B+ 


/o(980)A"+)B(/ (980) ->■ tt tt ) 


2.8 ± 0.6 ± 0.5 




2.8 ±0.8 


B(B+ 


-+ /o(980)iv + )B(/o(980) ->■ K+K~) 


9.4 ± 1.6 ± 2.8 


< 2.9 




B(B° 


-> fn(980)ftr°)B(fn(980) -> A"+K - ) 


7.0+55 

-3.0 




7.0+i^ 

— 3.0 


B(B+ 


— > /"nf980W+')B('/'(i(980 1 ) tt+tt - ~) 


< 1.5 




< 1.5 


B(B° 


— > /'n('980V)B('r'n(980') -> tt+tt -1 ) 


< 0.9 




< 0.9 


B(B° 


— i> fnOS0)ri)B! fn (980) — > tt+tt - ~) 


< 0.4 




< 0.4 


B(B + 


-+ a{j(980)K+)B(ao(980)° -> rpr°) 


< 2.5 




< 2.5 


B(B" 


— > a ( J(980)i4'+)S(ao(980) _ — > ^Tr - ) 


< 1.9 




< 1.9 


B(B + 


-> a^(980)A"°)B(ao(980)+ -> rpr+) 


< 3.9 




< 3.9 


B(B° 


-> ajJ(980)A"°)B(ao(980)° — > jjtt ) 


< 7.8 




< 7.8 


B(B+ 


-> ag(980)TT+)B(ao(980)° -> rpr°) 


< 5.8 




< 5.8 


B(B+ 


-> a+(980)Tr )B(a (980)° -+ rpr°) 


< 1.4 




< 1.4 


B(B° 


a+(980)Tr+)B(a (980)+ -> »7tt+") 


< 3.1 




< 3.1 


B(B° 


-> a+(1450)Tr+)B(a (1450)+ r?Tr+) 


< 2.3 




< 2.3 


B(B° 


-> a - (1450)if+)S(a (1450) - -> tttt - ) 


< 3.1 




< 3.1 


B(B+ 


-+ /o (1370)^+ W/ (1370) -> tt+tt - ) 


< 10.7 




< 10.7 


B(B+ 


/o(1370)tt+)B(/ (1370) -> tt+tt - ) 


2.9±0.5±0.5+°'g < 4.0 




< 4.0 


B(B+ 


-> / (1500)K+) 


17±4± 12 




17.0 ± 12.6 


B(B° 


-+ / (1500)^°) 


13.3+? ! ±3.2 

— 4.4 




13.3+^ 

— 5.4 


B(B+ 


-+ / (1710)J£'+)B(/o(1710) -> K+if - ) 


1.12 ±0.25 ±0.50 


1.7 ±1.0 


1.26 ±0.49 


B(B° 


-+ /o(1710)A' )B(/o(1710) ->■ X+K - ) 


4.4 ±0.7 ±0.5 




4.4 ±0.9 


B(B° 


-> /o(1710)K s )B(/o(1710) ->• X S X 5 ) 


0.50±°;g ± 0.11 




0.50+2-^ 


B(B+ 


-> i^*°(1430)Tr+) 


32.0 ± 1.2 +1 ° ® 

— 6.0 


51.6 ± 1.7±7g 


45.1 ± 6.3 


B(B° 


->■ i^p + (1430)Tr - ) 


29. 9+?/7 ± 3 - 6 " 


49.7 ± 3.8+g2 


qq c+3.9 

— 3.8 


B(B° 


-> ^:* (1430)tt ) 


7.0 ± 0.5 ± 1.1 6 




7.0 ± 1.2 


B(B+ 


->7fo°(1430)K+) 


< 2.2 




< 2.2 


B(B+ 


iC + (1430V) 


5.2 ± 1.9 ± 1.0 




5.2 ± 2.1 


B(B° 


K^ (1430)r7') 


6.3 ± 1.3 ± 0.9 




6.3 ± 1.6 


B(B+ 


-> ^+(1430)n) 


15.8 ± 2.2 ± 2.2 




15.8 ± 3.1 


B(B° 


-+ if*°(1430)r?) 

(J \ / // 


9.6 ± 1.4 ± 1.3 




9.6 ± 1.9 


B(B + 


— >• fn('980')/<"*+ , )BI' fnl'980 1 ) — i-Tr+Tr -, | 


42+06+03 




4.2 ± 0.7 


B(B° 


fn('980')A"* 01 lB(/'n(980') —> tt+tt -1 ) 


5.7 ± 0.6 ± 0.3 


< 2.2 


5.7 ± 0.7 


B(B+ 


— ¥ fn(98,0)o + )B( fnOSO) — > tt+tt - ") 


< 2.0 




< 2.0 


B{B° 


-> /'n('980')o ')B('/'n('980~) —5- TT+TT - ) 


< 0.40 


0.87 ± 0.27 ± 0.15 


0.87 ± 0.31 


B(B° 


— > /o(980)w)B(/o(980) tt+tt - ) 


< 1.5 




< 1.5 


B(B° 


—>■ fn('980~)ai')B('/'n('980~) — > tt+tt -1 ) 


< 0.38 




< 0.38 


B(B + 


— > i^Q + (1430)oj) 


24.0 ± 2.6 ± 4.4 




24.0 ± 5.1 


B(B° 


_> _ftT*°(1430)a;) 


16.0 ± 1.6 ± 3.0 




16.0 ± 3.4 


B(B+ 


-> i^*+(1430)cA) 


7.0 ± 1.3 ± 0.9 




7.0 ± 1.6 


B(B° 


-> ^*°(1430)<7i) 


3.9 ± 0.5 ± 0.6 




3.9 ± 0.8 


B(B° 


-s- K*°(1430)K*°) 




< 3.3 


< 3.3 


B(B° 


/sT* (1430)^* ) 


< 1.7 




< 1.7 


B(B° 


->iC*+(1430)p - ) 


28 ± 10 ± 5 ± 3 




28.0 ± 11.6 


B(B° 


-> ^*°(1430)p°) 


27 ± 4 ± 2 ± 3 




27.0 ±5.4 


B(B° 


-> /o(980)/o(980))B 2 (/ (980) ->■ tt+tt - ) 


< 0.19 


< 0.1 


< 0.1 


B(B° -> /o(980)if *°(1430))B(/ (980) tttt) 


2.7 ±0.7 ±0.5 ±0.3 




2.7 ±0.9 


B(B° 


->ii:* (1430)^ (1430)) 




< 4.7 


< 4.7 


B(B a -> K*°(1430)K* (1430)) 




< 8.4 


< 8.4 



a Thci£j is a new measurement of (27.8 ± 2 
b See 2|. 
c See 13 1 . 



5 ± 3.3) x 10~ 6 extracted from B° 
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TT by BaBar [2|. 



TABLE II: Experimental CP asymmetries (in units of %) of B decays to final states containing 
scalar mesons 111]. 



Mode 




BaBar 


Belle 


A 

Average 


r->4- 


• / (980)A + 


— 10.6 ± 5.0_^5 a 


i-r i-r \ r> r--|-4.6 

-7.7 ± 6.5I 2 .6 


n r+4 9 

-9- 5 -4.2 


B {) ->■ 


/ (980)A° 


-28 ± 24 ± 9 6 


-30 ± 29 ± 11 ± 9 b 


-29 ± 20 






-8 ± 19 ± 3 ± 4 c 


-6± 17±7±9 c 


-7 ± 14 


B + -) 


/ (l370)^+ 


72 ± 15 ± 14lg 




72 ± 22 


7-1 J- 


/o(1500)AT+ 


28 ± 26t 14 ' 




qq+30 
28 -29 


5+ -) 


ATq U (1430)7t + 


O.Z zt o.D 2 g 


7.6 ± 3.8^2:2 


r r+3.4 

0.0_2 2 




Ko + (1430)7T 


7 ± 14 ± 1 




7 ± 14 


7->n 


ATo u (1430)7r u 


— 15 ± 10 ± 4 




— 15 ± 11 


i? + -3 


ATq + (1430)77' 


6 ± 20 ± 2 




6 ± 20 




ATo u (1430)t7' 


-19 ± 17 ± 2 




-19 ± 17 


B + —? 


r * -4- / -1 1 fi n\ 

AT ^(1430)?7 


5 ± 13 ± 2 




5 ± 13 


£> u — > 


Aq(1430)?7 


6 ± 13 ± 2 




6 ± 13 


B + -; 


• fn(980 N )A'* + 


-15 ± 12± 3 




-15 ± 12 


S ^ 


/ (980)A*° 


7± 10 ± 2 




7 ±10 


5+^ 


(1430V 


-10±9±2 




-10±9 


5°^ 


AT *°(1430)w 


-7±9±2 




-7±9 


5+^ 


AT *+ (1430)0 


4± 15 ±4 




4 ±15 


5°^ 


A' O *°(143O)0 


20 ± 14 ± 6 




20 ± 15 



"This data is from the measurement of B + K + n + ir . BaBar's measurements of B + — > K + K + K 
and B+ -» A+7r 7r° yield A CP (B+ -> / (980)A+) = -(8±8±4)% 0] and A CP {B+ -> / (980)A+) = 
(18 ± 18 ± 4)% 5], respectively. 

fc From B" -> K+R-Rs- 

c From B° 7r + 7i-2fs- 



we assume in scenario 2 that Aq (1430), ao(1450), /o(1500) are the lowest lying qq resonances and 
the corresponding first excited states lie between (2.0 ~ 2.3) GeV. Scenario 2 corresponds to the 
case that light scalar mesons are four-quark bound states, while all scalar mesons are made of two 
quarks in scenario 1. Phenomenological studies in 0,0] imply that scenario 2 is preferable, which 
will be also reinforced in this work. Indeed, lattice calculations have confirmed that ao(1450) and 
A"q(1430) are lowest-lying P-wave qq mesons |2{j], and indicated that /o(500) (or a) and Aq(800) 



(or k) are S-w&ve tetraquark mesonia 30]. 1 



However, a recent lattice calculation 



31] leads to an opposite conclusion. 
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TABLE III: The scalar decay constant fs (in units of MeV) and Gegenbauer moments B%, B% and 
in scenario 1 (left) and scenario 2 (right) at the scale [i = 1 GeV obtained using the QCD sum 
rule method [6j]. Decay constants and Gegenbauer moments for excited states in scenario 2 are not 
listed here. 





Is 


B x 


B 3 


fs 


Bi 


B 3 


/o(980) 


370 ± 20 


-0.78 ± 0.08 


0.02 ±0.07 








a (980) 


365 ± 20 


-0.93 ±0.10 


0.14 ±0.08 








Kl (800) 


340 ± 20 


-0.92 ±0.11 


0.15 ±0.09 








/ (1500) 


-255 ± 30 


0.80 ± 0.40 


-1.32 ±0.14 


490 ± 50 


-0.48 ±0.11 


-0.37 ±0.20 


a (1450) 


-280 ± 30 


0.89 ±0.20 


-1.38 ±0.18 


460 ± 50 


-0.58 ±0.12 


-0.49 ±0.15 


if *(1430) 


-300 ± 30 


0.58 ±0.07 


-1.20 ±0.08 


445 ± 50 


-0.57 ±0.13 


-0.42 ±0.22 



TABLE IV: Form factors of B 
front model 



ao ( 1450), Kq{ 1430) transitions obtained in the covariant light- 



321 ] in scenario 1 (upper entry) and scenario 2 (lower entry). 



F F(0) a b 


F F(0) a b 


^(1450) 0.26 1.57 0.70 
0.21 1.66 1.00 

F BK*(U30) 1 

0.21 1.59 0.91 


^0(1450) Q26 

0.21 0.73 0.09 

F BK*(1430) 

0.21 0.59 0.09 



A. Decay constants and form factors 

Decay constants of scalar, pseudoscalar and vector mesons are defined as 

(S(p)\q2l»qi\ty = fsP», (S\ qiq\ |0) = m s fs, 
(P(p)\q2j^qi\0) = -ifpPv, (^(p)|<?27/Wi|0) = forays*^ 
(V(p,e*)\q2^ qi \0) = f^sl -Pvsl). (2.1) 

For scalar mesons, the vector decay constant fs and the scale-dependent scalar decay constant fs 
are related by equations of motion 

vsfs = fs, with [i s = — — rr> ( 2 - 2 ) 

m 2 W - mi(/i) 

where and mi are the running current quark masses and mj is the scalar meson mass. For 
the neutral scalar mesons /o, a[j and a, fs vanishes owing to charge conjugation invariance or 
conservation of vector current, but the quantity fs = fJ-sfs remains finite. It is straightforward to 
show from Eq. (|2.ip that the decay constants of the scalar meson and its antiparticle are related 
by 

fs = fs, fs = ~fs- (2-3) 
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Indeed, from Eq. (|2.2p we have, for example, 

m d (fi) - m u (/i) 

, ! /' ) — J,,., 



m u {n) - m d {n) 



ao m ao 

Therefore, the vector decay constants of ag and ag" are of opposite sign. 

In we have applied the QCD sum rule method to estimate the decay constant fg for various 
scalar mesons as summarized in Table IIIH Note that a recent sum rule calculation in [25] yields a 
smaller fg in scenario 2 for S = /n (1710), an (1450) and ETq(1430). In this work we shall use the 



(2.4) 



values of fy and fy taken from 33J. For the decay constants and /* (/) of the rj and rf mesons 
defined by 



(2.5) 



we shall follow the results of 
For the B — >• P and B — > 



34] 



V transition form factors defined in the conventional way 35j, we will 
use the results obtained using the QCD sum rule method [36]. Form factors for B — > S transitions 
are defined by [32| 



(S(p')\AjB(p)) 



where P M = {p + p')^ and q^ = (p 



P„ 



m 



B 



mt 



% ) F? S tf) + 



mi 



mt 



(2.6) 



parameterized in a 3-parameter form 

F(q 2 ) - 



p')a- The momentum dependence of the form factor is usually 
F(0) 



(2.7) 



1 - a(q 2 /m 2 B ) + b(q 4 /m%) ' 
The parameters F(0), a and b for B — ^ S transitions are summarized in Table ITVl obtained using 
the covariant light-front quark model [32|. Form factors are also available in other approaches, such 
as light-cone sum rule 043 and pQCD QQ- 

In general, form factors calculated by sum rule 
and pQCD methods are larger than that obtained using the quark model. For example, (0) 
is of order 0.26 in the covariant light-front quark model 32] , while it is found to be 0.45 25], 0.49 
221 . |2J] in the sum rule method and 0.60 [261 ] and 0.76 iJ] in pQCD (all evaluated in scenario 2). 
We will come to this point later. 



B. Distribution amplitudes 



In general, the twist- 2 light-cone distribution amplitude (LCDA) of the scalar meson $5 has 
the form 



$s(x,(J>) = /sr6x(l - x) 



1 + ^^^)^(2,-1) 



771=1 



(2.8) 



3/2 

where B m are Gegenbauer moments and Cm are Gegenbauer polynomials. The general twist-3 
LCD As are given by 



*sW = Is 



1+ ]T a m (fi)CU 2 {2u-l) 



777=1 



®g(x) = fg6x(l-x) 



l+Y j bMC l J l 2 {2u-l) 



777=1 



(2.9) 
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Since us = 1/Bq S> 1 and even Gegenbauer coefficients B m are suppressed, it is clear that the 
twist-2 LCDA of the scalar meson is dominated by the odd Gegenabuer moments. In contrast, 
the odd Gegenbauer moments vanish for the n and p mesons. The Gegenbauer moments B\ and 
B3 in scenarios 1 and 2 obtained using the QCD sum rule method 0] are listed in Table IIII[ The 



Gegenbauer moments 0124 and 61,2,4 for twist-3 LCDAs have been computed in 25l. l37|. 

7° 

'0 



Since the decay constants vanish for the neutral scalar mesons /o, aj] an d it follows from Eq. 



(g3D that 

<D 5 (x, u) = f s 6x(l ~x)J2 B m (u) Cli 2 {2x - 1) (2.10) 

m=l 

for these neutral scalar mesons. 

As stressed in [g], it is most suitable to define the LCDAs of scalar mesons including decay 
constants. However, it is more convenient in practical calculations to factor out the decay constants 
in the LCDAs and put them back in the appropriate places. In the ensuing discussions, we will use 
the LCDAs with the decay constants fsifsifvi fyifp being factored out. 

C. Mixing angle between /o(980) and /o(500) and between 77 and rj 

In the naive 2-quark model with ideal mixing for /o(980) and /o(500), /o(980) is purely an ss 
state, while /o(500) is a nn state with nn = (uu + dd)j\[2. However, there also exist some experi- 
mental evidences indicating that /o(980) is not purely an ss state. For example, the observation of 



Y(J/ip — > foui) ^r(J/i[) — y fo4>) 2a] clearly shows the existence of the non-strange and strange 



quark content in / (980). Therefore, isoscalars /o(500) and /o(980) must have a mixing 

( l/o(980)) \( cos 9 sm6 \ ( \nn) \ 
Vl/o(500)) J \-sm6 cosO) \ \ss) J ' 

Various mixing angle measurements have been discussed in the literature and summarized in (o. 39 1. 
A recent measurement of the upper limit on the branching fraction product B(B° — > J/^/o(980)) x 



/3(/ (980) -> vr+TT-) by LHCb leads to \9\ < 30° [4C]. 



For the i] and rj' mesons, it is more convenient to consider the flavor states qq = (uu + dd)/y/2, 
ss and cc labeled by the rj q , rj s and r/J?, respectively. Neglecting the small mixing with ryj?, we write 

\v} \ _f cos -sin0\ / \r) q }\ 
\ sin^ cos0 I \ \rj s ) J 



(2.12) 



where cj) = (39.3 ± 1.0)° 41[ is the r] — rj' mixing angle in the r\ q and ij s flavor basis. 



III. DECAY AMPLITUDES IN QCD FACTORIZATION 



We shall use the QCD factorization approach 27|, |28| to study the short-distance contributions 



to the B -> SP,SV decays with S = f (980), a (980), a ( 1450), ^(1430). In QCD factorization, 
the factorizable amplitudes of above-mentioned decays can be found in [y] and 0]. However, the 
expressions for the decay amplitudes of B — > JqK, a^n, QqK involving a neutral / or oq given in [6] 
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are corrected in Appendix A as some factorizable contributions were missed before. The effective 
parameters a% with p = u,c appearing in Eq. (|Aip can be calculated in the QCD factorization 
approach (27]]. In general, they have the expressions 

a^M 1 M 2 )=^c l + C ^N l (M 2 ) + C ^^[v t (M 2 ) + ^H l (M 1 M 2 )} +if(M 2 ), (3.1) 

where i = 1, • • • , 10, the upper (lower) signs apply when i is odd (even), a are the Wilson coefficients, 
C F = (A 2 - l)/(2N c ) with N c = 3, M 2 is the emitted meson and Mi shares the same spectator 
quark with the B meson. The quantities Vi{M 2 ) account for vertex corrections, Hi(M\M 2 ) for hard 
spectator interactions with a hard gluon exchange between the emitted meson and the spectator 
quark of the B meson and Pi{M 2 ) for penguin contractions. The expression of the quantities 
Ni{M 2 ) reads 

, . f 0, i = 6, 8 and M 2 = V, . . 

N i {M 2 ) = \: ' (3.2) 

I 1, else. 

The explicit expressions of Vi(M), Hi(MiM 2 ), and weak annihilation contributions described by 
the terms bi and 6i,E\v are given in 0] and |t| for B — > SP and B — > SV, respectively. 2 

Power corrections in QCDF always involve troublesome endpoint divergences. We shall follow 



27] to model the endpoint divergence X = Jq 1 dx/(l — x) in the annihilation and hard spectator 
scattering diagrams as 

X A = In (1 + p A e«*), X H = In (1 + pse**), (3.3) 

with Ah being a typical scale of order 500 MeV, and pa,h, 4>a,h being the unknown real parameters. 

In principle, physics should be independent of the choice of p, but in practice there exists some 
residual p dependence in the truncated calculations. However, we found that sometimes even the 
decay rates without annihilation are sensitive to the choice of p. For example, we found that the 
measured branching fractions of B — > Kq (1430) (77, 7/) cannot be accommodated for p = 
Indeed, this observation also occurs in our previous study of B — > VV decays [42]. We found 
that if the renormalization scale is chosen to be p = m\,{mb)/2 = 2.1 GeV, we cannot fit the 
branching fractions and polarization fractions simultaneously for both B — > K*(p and B — > K* p 
decays. Therefore, we will confine ourselves to the renormalization scale p = mb(mb) in the ensuing 
study. Note that the hard spectator and annihilation contributions should be evaluated at the hard- 



collinear scale ph = \fpKh with A/j « 500 MeV [27]. 



As discussed in [6|] and [7J, scenario 2 in which the scalar mesons above 1 GeV are lowest lying qq 
scalar state and the light scalar mesons are four-quark states is preferable, while all scalar mesons 
are made of qq quarks in scenario 1. It is widely believed that the /o(980) and the ao(980) are 
predominately four-quark states, but in practice it is difficult to make quantitative predictions on 
hadronic B — > SP, SV decays based on the four-quark picture for light scalar mesons. Hence, we 
shall assume scenario 1 for the /o(980) and the ao(980) in order to apply QCDF. 



In Eq. (4.8) of [6], the second term of the annihilation amplitude A3 for M\Mi = SP should have an 
identical expression, including the sign, as that for M\Mi = PS. Note that the expression was correct in 



the original archive version, arXiv:hep-ph/0508104 
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TABLE V: Branching fractions (in units of 10 -6 ) of B decays to a scalar meson and a pseudoscalar 
meson. The predicted rates of B — > /o (980) if, /o(980)7r are for the /o(980) — /o(500) mixing angle 
9 = 17°. We work in scenario 1 for the light scalar mesons /o(980) and ao(980) and scenario 2 
for the scalar mesons ao(1450) and if q (1430); see the main text for explanation. Experimental 
results are taken from Table HI We have used B(/o(980) — > 7T + 7r~) = 0.50 and ,6(ao(980) — > ttt]) = 
0.845 db 0.017 to obtain the experimental branching fractions for /o(980)P and ao(980)P. For 
comparison, predictions based on the pQCD approach are also exhibited. 



Mode 




QCDF (this work) 


pQCD [11-14, 11] 


Expt 


B~ - 


/o(980)Jf- 


1 R 1 +1-9+1-2+30.8 
1U - J --1. 8-1. 8-11.0 


16- 18 a 


18.8^ 


j5° -» 


io(980)if° 


i a q + 1. 7+1. 1+28.6 
L *'°—l. 6—1. 6-10.2 


13- 16 a 


14.0 ± 1.8 


B~ - 


/o(980)vr- 


n oR+0. 04+0. 05+0. 18 
u - zo -0. 03-0. 04-0. 12 


2.5 ± 1.0 


< 3.0 


3° — ► 


/o(980)^° 


n ns+ - 01 +°- cll + - 08 

u - u °-0. 01-0. 01-0. 03 


0.26 ±0.06 




B~ - 


a%(980)K- 


n OA+0- 08+0. 45+1. 02 
06-0. 13-0. 07 


q r+0. 4+0. 4+1.0 
°- -0. 4-0. 6-1.0 


< 3.0 


B~ - 


a^{980)K° 


n ns+°- 12 +°- 58 + 212 

U.UO_q q7_q .08-0 .02 


fi q+0. 8+1. 1+2.0 
"•''—0.7—1.1—1.7 


< 4.6 


B° -> 


a£(980)K- 


n OA+0. 05+0. 62+2. 27 
u - oq: -0. 05-0. 13-0. 00 


Q 7+I.I + I.6+2.7 
J - ' -1.0-1.4-2.2 


< 2.2 


B° -> 


a° (980)K° 


n nc+0. 05+0. 28+0. 86 
u - uo -0. 03-0. 04-0. 00 


1 7+0.5+0.7+1.1 
' -0.5-0.8-1.1 


< 9.2 


B~ - 


a[j(980)7r- 


a Q+0. 4+1. 2+0.5 
J -0.3-l.l-0.7 


9 o+O.O+O.O+O.O 
z -°-0.8-0.9-0.6 


< 6.9 


B~ - 


a (980)tt° 


n 7n +0. 22+0. 09+0. 22 
u - ' u -0. 16-0.08-0. 14 


n /1 +0.00+0.00+0.00 
u -* 1 -0.13-0. 14-0. 12 


< 1.7 


~B° — ► 


a([(980)7r- 


K o+0. 3+1. 4+1.0 
J -°-0.4-1.3-0.4 


n oe+0. 10+0. 14+0.01 
u -°°-0. 09-0. 14-0. 00 




B -» 


ag (980)tt+ 


n co+0. 11+0. 09+0. 63 
U.OO_q og-0 .08-0 .22 


n ri +0.05+0.09+0.07 
U- 01 -0.O6-0. 09-0. 06 




B° -» 


a^(980)7rT 


z, y+0.6+1.3+1.8 
*• ' -0.6-1.1-0.5 


n QQ+ .l +°15+0 02 
u - yo -0. 10-0. 14-0.00 


< 3.7 


B -» 


ag(980)7T° 


i n+0.5+0.1+0.2 
-0.3—0.1 — 0.1 


n ri +0.08+0.09+0.00 
u,d — 0.07— 0.09— 0.00 




s- - 


a°,(1450)if- 


9 9+0.7+2.3+7.7 
0.5— 1.0— 1.9 






B" - 


0(7(1450)^° 


/I 9 + 1.6+4.9+18.1 
*- z -1.2-2.1-4.2 






j5° -» 


a+(1450)if" 


q c + 1. 0+4.5+16. 9 
°' J -0.8-1.9- 3.3 




< 4.7 


B° -> aSj(1450)i?° 


i n+0.8+2. 1+7.5 
1- — 0.6-0.9— 1.9 






s- - 


ag(1450)7r- 


K -,+0.5+1.2+1.3 
tJ " L — 0.4— 1.1— 1.3 






s- - 


0(7(1450)77° 


9 -, +0.7+0.2+0.8 
z - x -0. 5-0.2-0.6 






B°-> 


fl([(1450)7r- 


9 r+0.5+0.9+3.7 
z - d -0.5-0.8-0.3 






j5° -> 


flg (1450)7T + 


n 7/ i+0.20+0.19+2.92 
u - '^-0.16-0.17-0.51 






j5° -» 


aJ(1450)7rT 


i q+0. 7+0. 7+8. 3 
i '°-0. 6-0. 6-0.1 




< 3.5 


B° — >■ 


ag(1450)7r° 


o o + l. 4+0. 3+2. 8 
°-°-1.0-0.3-1.3 






s- - 


i?,2; (1430)7r- 


19 q+4.6+4. 1+38.5 
lz - y -3. 7-3. 4-9.1 


30.91 1 !;! b 


45.1 ±6.3 


s- - 


K*-{U30)n° 


7 A +2. 4+2. 1+20.1 
' '^-1.9-1.8- 5.0 


91 R+8.5 6 
zl '°-6.6 




i?° — > 


K*-(U30)tt+ 


1 •} c+4.5+4.1+38.3 
i0 -°-3.6-3.5- 9.5 


SI.et- 4 » 


qq c+3.9 


i?° — > 


if *°(1430)7r° 


r R +2.6+2.4+18.8 
°-°-1.3-1.2- 3.9 




7.0 ± 1.2 


B- - 


iC*- (1430)7/ 


i 7 Q+3.9+8.3+ 9.1 
1 ' - y -3.4-5.3-12.3 


qq p + 13. 5 + 1. 1 + 7. 7 + 8. 2 

9.0-1.1-7.0-7.3 


15.8 ±3.1 


B- - 


K*-(U30) V > 


Q q+4.7+4.0+51.6 
y -°-3.6-4.4- 8.0 


77 r + 15. 8+6. 2+21. 0+18.0 
' ' -°-10. 8-5. 8-16. 5-16.1 


5.2 ±2.1 


B° — > 


i? *°(1430)?7 


10 1+3.6+7.6+ 9.1 
iu ' 1 -3.1-4.9-11.7 


90 4+11.6+1.4+6.4+6.9 
zo -^- 7.8-1.4-5.9-6.2 


9.6 ± 1.9 


i?° -> 


if *°(1430)r/ 


7+4.4+3.7+48.7 
°-' -3.3-4.1- 7.5 


7A 9+15.0+6.4+20.5+17.2 
'^• z - 10.3-5.7-16.2-15.5 


6.3 ± 1.6 



11]. 



Tor the mixing angle 140° < 6 < 165° 
^Results based on the new Gegenbauer moments9obtained in [14j . 
B -> isT^(1430)7r, see \P 



For previous pQCD calculations of 



TABLE VI: Same as Table IVl except for B decays to a scalar meson and a vector meson. 



Mode 



QCDF (this work) 



pQCD [15,16, 18- 



Expt 



B° 
B° 
B~ 
B~ 
B° 
B° 
B° 
B~ 
B° 



B° 
B° 
B~ 



B° 
B~ 
B~ 
B° 
B° 
B~ 
B° 
B~ 
B° 



*o 



B~ -> / (980)ir* 
/o(980)A 
> /o(980)p- 
/o(980)p° 
/o (980)^ 
• a°,(980)A: 
B- -»■ a (980)A 



B° 
B- 
B° 
B° 
B- 



*o 



ag(980)^*° 
• ag(980)p- 

► (980)p° 
0^(980)^- 

(980)p+ 
a , (980)/ 

► (980)w 
ag(980)a; 

► ag(1450)A'*- 
B- -> a (1450)A 

^(1450)A*- 



"F7*0 



ag(1450)A 
► ag(1450)p" 
B- -> ao(1450)p° 
B° 
B° 
B° 



a^(1450)p- 
Oq (1450)/9 + 
a[](1450) / 9 



B- -»■ (1450)w 
a°,(1450)u; 

► #S°(1430)/T 

► A*-(1430)p° 
A*-(1430)p+ 
A > ^°(1430)p° 

► A*"(1430)w 



A *°(1430)w 
► tf -(1430)< 
A^°(1430)</) 



1n 7+0.3+0.2+4.8 
iu - ' -0.3-0.2-0.5 



9.1 



+ 1.0+1.0+5.3 
-0.4-0.5-0.7 



n oc+0.01+0.08+0.09 
u - od -0.01-0.07-0.08 

n no+O-OO+o.00+0.03 
u - uz -o.oo-o.oo-o.oo 
n no +o.oo+o.oo+o.03 

u - uz -0.00-0.00-0.01 
1 +0.3+Q.4+2.i 



2.4: 

8.4 



-0.3-0.4-0.0 
+1.3+1.2+8.2 



1.1-1.1-0.2 
r fi +l. 2+1. 4+7.9 
J -°-0.2-0.4-0.6 
2.5+2.7+2.8 



3.2 
6.8 
2.1 
22.5 
0.60 
1.3 
1.0 



2.9-2.7-0.0 
+0.6+1.7+0.9 
0.5-1.5-0.3 
+0.6+0.0+0.9 
0.5-0.0-0.2 
-4.0+6.8+5.0 
-0.1-2.7-2.3 
+0.14+0.05+0.47 
0.03-0.05-0.07 
+0.8+0.1+0.4 
0.4-0.1-0.0 
+0.4+0.0+0.5 
0.3-0.0-0.2 



n ,.+0.3+0. 0+0.2 
u -*-0.1-0. 0-0.0 
9 c;+0.3+0.4+5.6 
z - d -0.3-0. 4-0.1 
o „+l. 4+1. 3+14.6 
°-^-l. 1-1.2- 0.4 
o q+1.9+3.0+3.9 
°- y -0.3-l. 4-0.5 
o ,,+0.8+0.9+6.7 
°-*-0.2-0. 3-0.3 

7+0.6+1.4+1.9 
°' -0.5-1.1-0.3 

+0.8+0.2+2.3 
0.6-0.2-0.4 
2.4+4.5+11.2 
0.4-2.3- 1.0 
+0.3+0.1+1.3 
0.1-0.1-0.1 
+2.7+0.5+1.1 
1.3-0.2-1.2 

1 r+0.6+0.0+1.5 
1 - d -0.4-0. 0-0.3 

2+1.1+0.2+1.7 
'-0.5-0.1-0.3 
+34.5+35.2+41.6 
-35.8-35.3-51.0 
+3.7+0.4+6.7 



3.2 
11.2 
1.2 
4.4 



1.. 
39.0 



14-8_3 .2-0.4-0.2 

36 3+ 8 - 5 +°- 8 + 16 - 7 



23.4" 



-7.4-0.8- 0.1 
-5.1+0.6+9.4 
-4.5-0.5-0.4 
oi c+5. 8+0.5+9.8 
zi -°-4.9-0.5-3.4 
01 q+5.9+0.6+10.5 
zi - y -5. 0-0.6- 3.3 
q O+0.7+0.1+11.2 
°-°-0. 6-0.1- 1.6 
q 7+O.8+O.I+5.6 
°-' -0.6-0.1-0.1 



11.7 ~ 14.6 
11.2 ~ 13.7 

n 7+0.1+0.2+0.2 
u -' -0.0-0.1-0.1 
n qq+0.04+0.07+0.06 
u -°°-0.03-0.05-0.06 



0.34 



+0.03+0.06+0.06 
-0.04-0.06-0.05 



7 n +0. 9+1.6+1. 7+0.2 
' - u -0. 7-1. 1-1. 4-0.0 
q ,-,+0.2+0.2+0.7+1.2 
°- u -0.1-0.1-0.6-0.7 
o o+0. 3+0.1+0. 7+0.8 
z -°-0.3-0.0-0. 5-0.6 
'.1+0.0+0.3+0.5 



1 4+0-U 



0.1-0.3-0.4 
7Q q+33.2+18.2+5.5+3.9 
<y-d-22. 1-16. 3-4.7-3. 5 
+0.9+0.4+0.1+0.4 



1.9 
184.3 
3.6 
7.2 



0.5-0.4-0.0-0.3 
+69.3+42.2+13.0+13.7 
47.9-37.9-13.3-14.1 
+0.7+0.8+0.8+0.4 
0.6-0.8-0.8-0.3 
+2.1+1.7+1.1+1.6 
1.6-1.3-0.9-1.3 
q 3+0.4+0.1+0.1+0.1 



2.1 
12.1 
8.4 



0.1-0.0-0.0-0.0 
+0.7+0.4+0.2+0.1 
0.6-0.5-0.2-0.2 
+2.8+3.9+0.5 
0.0-3.1-0.5 
+2.3+3.3+0.9 
0.0-3.2-0.7 



qn c+2. 7+3. 5+0.3 
lu -" J -0.0-2. 6-0.3 



4.8; 



-1.1+1.0+0.3 



-0.0-1.0-0.3 
7 .+2.1+3.0+0.9 
' -^-1.5-2.3-0.4 
3+2.1+3.6+1.2 
-2.0-2.9-1.0 
oc fi +6.2+0.9+12.1 
zd - u -5.4-0.8- 6.5 
q fi +5.6+0.8+10.9 
ZO - U -5.0-0.6- 5.8 



9.3] 



11.4 ±1.4 
8.4 ± 1.4 

< 4.0 
1.7 ±0.6 

< 3.0 



28.0 ± 11.6 
27.0 ±5.4 
24.0 ±5.1 
16.0 ±3.4 
7.0 ± 1.6 
3.9 ±0.8 



10 



TABLE VII: Branching fractions (in units of 10~ 6 ) of B — > SP (upper) and B — > SV (lower) decays 
with S = ao(1450), Kq (1430) in QCD factorization. Experimental results are taken from Table HI 
The scalar mesons ao(1450) and .Kg (1430) are treated as the first excited states of low lying light 
qq scalars ao(980) and .K"q(800), corresponding to scenario 1 as explained in the main text. 



Mode 



Theory Expt 



Mode 



Theory Expt 



w 
if 
if 
if 
if 
if 
if 

if 
if 
if 



B- 
B' 
B- 
B~ 



B' 
B' 
B' 
B~ 
B~ 
B~ 
B- 
B' 
B' 



a%(U50)K- 
(Iq (1450)A° 
ag(1450)vr~ 
(1450)tt° 



^°(1430)^- 
A *-(1430)tt 
K *-(1430)ry 
A'*- (1430)7/ 



0.7 
2.7 



+0.7+0.3+3.1 



-0.4- 
+2.8- 



0.1-0.5 
1.2+9.6 



-1.7-0.6-2.7 
9 7 +0.1+0.6+0.4 
z - ' -0.1-0.5-0.5 
a+0. 2+0.0+0. 2 



0.4' 



0.2-0.0-0.2 



1.3 
0.3 
5.4 

6.2 



+0.7+1.0+15.9 
0.5-0.7- 1.1 
+0.2+0.3+6.4 
0.1-0.2-0.3 
+1.9+1.2+3.7 
1.6-1.0-4.4 
+ 1.0+1.0+35.8 
-0.8-1.0- 6.1 



45.1 ±6.3 

15.8 ±3.1 
5.2 ±2.1 



W 

if 

if 

if 

if 

if 

2° 

if 

if 

if 



a$(1450)K- 

a° (U50)K° 

0^(1450)^- 

ao(1450)vr+ 

a^(1450)7r =F 

ag(1450)rr 

A *~(1430)^ + 

Ao°(1430)7T 
IF* 

K o V 



1.9 
0.9 



+ 1.9+0.9+7.8 
1.1-0.4-1.4 
+1.0+0.4+3.6 
0.6-0.2-0.9 



< 4.7 



9+2.0+1.8+4.5 
ii - z -l. 6-1. 7-5.2 



0.02 
11.9 

1.3 
1.6 
1.3 

5.8 
5.9 



+0.05+0.01+0.75 
0.00-0.00-0.01 
+2.7+1.8+5.6 
2.3-1.6-3.7 
+0.7+0.1+2.0 
0.5-0.1-1.0 
+0.8+1.0+15.7 
0.6-0.8- 1.3 
+0.7+0.8+9.4 
0.3-0.5-1.3 
+2.1+1.4+4.9 
-1.7-1.1-5.2 
0.9+0.9+33.6 
-0.8-0.9-5.8 



< 3.5 

32.1 ±3.7 
7.0 ± 1.2 
9.6 ± 1.9 
6.3 ± 1.6 



B- 
B- 
B~ 



*o 



1.3 



ag(1450)A'*- 0.7; 

(1450) A 
a[j(1450)p~ 
ao(1450)p° 
ao(1450)a; 



-0.1+0.2+1.7 
-0.1-0.2-0.0 
+0.2+0.5+4.2 
-0.1-0.4-0.2 
g g+1. 3+1. 8+0.2 



-1.1-1.' 



-1.1 



0.2 
0.5 



+0.1+0.0+0.6 
-0.1-0.0-0.1 
+0.1+0.1+0.4 
-0.1-0.1-0.1 



4.81°; 
o.6±°; 

+4. 
-1. 
+0. 



iC *~(1430)</> 2.8 
A*°(1430)p- 
A *"(1430)p° 
X~o~(1430)a; 



+0.5+0.0+5.0 
0.4-0.0-0.5 
i o 9+ 9.9+10.4+10.8 
io - z -10.8-10.4-15.1 
o A +1.7+0.2+2.1 
°-^-l. 5-0.2-0.1 
in c+3.0+0.2+3.8 
J-U-O_2.4-0.2-L4 



7.0 ± 1.6 



24.0 ±5.1 



aj(1450)__*- 
a° (U50)K*° 
a£(U50)p- 
a (U50)p + 
a%(U50)p° 
ag(1450)w 
Kf (1430)0 
K*~(U30)p + 10.9 



17.3 
0.2 



o-i±S 



2.5 



Ao°(1430)p° 



4.11} 
9.3i' 



9+2.1+7.3 
1-1.1-0.4 
1+0.3+1.9 
0-0.2-0.1 
2+5.6+2.3 

0- 2.4-6.0 
1+0.0+0.3 

1- 0.0-0.0 
6+0.2+0.3 
8-0.1-0.8 
1+0.0+0.5 
0-0.0-0.0 
5+0.0+4.7 
3-0.0-0.5 
5+0.3+5.5 

2- 0.3-0.0 
1+0.2+2.6 
0-0.2-0.1 
7+0.3+3.9 
2-0.3-1.3 



3.9 ±0.8 

28.0 ± 11.6 
27.0 ±5.4 
16.0 ±3.4 



IV. RESULTS AND DISCUSSIONS 

The calculated branching fractions and CP asymmetries of B — > SP and B — > SV decays are 
summarized in Tables IVHlXl The defaulted values of the parameters pa,h and (j>A,H introduced in 
Eq. (|3.3p are set to zero; that is, the central values (or "default" results) correspond to pa,h = 
and 4>a,h = 0. The first theoretical error shown in the Tables for QCDF results is due to the 
variation of B\^ and fs, the second error comes from the uncertainties of form factors and the 
strange quark mass, while the third error from the power corrections due to weak annihilation and 
hard spectator interactions. 

In order to compare theory with experiment for decays involving /o(980) or ao(980) and ao(1450), 
we need an input for £>(/o(980) — > 7r + 7r~) or B(oq — > nrj). To do this for /q(980), we shall use the 
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TABLE VIII: CP asymmetries (in units of %) of B decays to a scalar meson and a pseudoscalar 
meson. We work in scenario 1 for the light scalar mesons _/o(980) and ao(980) and scenario 2 for 
the scalar mesons ao(1450) and i£g(1450). Experimental results are taken from Table HH 



Mode 




Theory 


Expt 


Mode 




Theory 


Expt 


B~ -> 


/o(980)X- 


i x+0.2+0.5+1.1 
J-*— 0.2-0.4— 1.4 






/o(980)if 


9 q+0.4+0.5+4.8 
z - y -0.3-0.2-2.6 


-14± 12 a 


B~ -> 


/o(980)tt- 


in n+3. 1+2. 0+36.1 
lw - u -2.6-2.4-28.4 




n 

B°^ 


/o(980)tt 


qi n+3.5+3.1+65.5 
J4: - u -3.9-2.9-47.2 




B~ -> 


ag(980)if- 


i o o+2.7+7.3+56.8 
"■^ — 2.7—6.1—40.1 




B° -> 


a£(980)K- 


i q 1 +5.3+10.1+75.2 
J - a -- L -6.4- 8.6-41.2 




B~ -> 


ap (980)X 


n cq+0. 74+54. 22+2. 49 
U.DO_ 57 _ 0.49-2.86 




n 

B°^ 


a°(980)if 


7 fi + 1-9+58.2+35.6 
'•°— 3.9- 3.0-30.5 




B- -+ 


a[j(980)7r- 


o n+0.6+0.0+0.9 
°-°-0.6-0.0-0.8 




B°^ 


a+(980)7r~ 


i -,+0.1+0.1 + 16.5 
•— L -0.1-0.0-16.9 




B~ -> 


a (980)7T° 


qo O+5.6+3.2+20.3 
JZ -°-6.2-2.8-21.6 




B°^ 
B° -> 


a (980)tt+ 
a[j(980)7r° 


r n + l. 0+0.6+16. 6 
d - u -0.9-0.5-15.5 
qn n+4.6+1.9+30.5 
ou ' y -5.6-1.8-24.7 




B~ -> 


ag(1450)iir- 


1 7Q+0. 69+0. 64+16. 61 
i - ' J -0. 82-0. 06-17. 75 




B° -> 


a+(1450)JsT- 


n fic+0.41+0. 68+11. 29 
u - UJ -0. 54-0. 23-23. 80 




B~ -> 


n 

Oq (1450)X 


n 99+0.04+0.06+ 0.59 
u - zz -0. 03-0. 03-47. 94 




n 

B°^ 


a°(1450)7T 


n 7A +0.32+0.23+ 8.2 
u - '^-0.39-0.05-20.7 




B~ -+ 


a[j(1450)7r~ 


^ r. ^ -1-0 9+0 04-4 9 
— lO.±_ g_ o_ 4 3 




-ft 

B -+ 


a + (1450)7T" 


^ +0 T74-0 18+8? 8 
— u - zl -0.46-0.28-83.6 




B" -> 


a (1450)tt 


9/1 fi+3. 4+2.4+25. 3 
z ^-°-3. 8-2. 1-23.5 




B°^ 
B° -> 


a "(1450)7r+ 
a°,(1450)77° 


r 9+3.4+2.3+62.0 
°- z -2.9-1.7-66.9 
i q o+3. 1+2. 1+35.7 
10 -°-3. 8-1. 9-20.0 




B" -> 


Fo°(1430)7r- 


i q+0. 1+0.0+5.9 
1,o -0. 1-0.0-4.8 


5.5^2'2 


B°^ 


X *-(1430)7r+ 


n 91 +0.06+0.05+7.05 
u ' zl -0. 06-0. 04-7. 21 


7± 14 


B" -+ 


(1430)tt° 


q n+0.4+0.5+10.7 
°- u -0.4-0.4- 7.8 




B°^ 


F; (1430)7r 


1 q+0. 4+0. 4+12.0 
i - J -0.5-0.5-14.9 


-15 ±11 


B" -+ 


^*-(1430)r? 


i (-+0.1+0.2+X.7 
•'••"-0.1-0.2-1.9 


-19 ±17 


B°^ 


i?*°(1430)ry 


9 n +0. 1+0.2+2.3 
z - u -0. 1-0.2-1.3 


6 ±13 


B~ -+ 


K*- (1430)77' 


i 9+0.3+0.1+8.7 
••••^-0.2-0.3-6.3 


6 ±20 


B° -> 


K* ° (1430)7?' 


1 9+0.1+0.1+17.9 
i - z -0. 1-0.4-13. 5 


-19 ±17 



"This is the naive (uncorrelated) average of the direct CP asymmetries —0.29 ± 0.20 obtained from B° 
K+K-Ks and -0.07 ± 0.14 from B° -> 77+77-^5 (see TablcHJ). 



BES measurement 43f] 



r(/ (980) -> 7777) 



= 0.75; 



-0.11 
-0.13 • 



r(/ (980) 7777) +r(/ (980) -> iYiY) ' 1 U) 

Assuming the dominance of the /o(980) width by 7777 and KK and applying isospin relation, we 
obtain 3 

B(/ o (980) -> 77+77-) = 0.501°$ , B(/ o (980) -> if+tf-) = 0.125l°;°22 • ( 4 - 2 ) 

For ao(980), we shall apply the Particle Data Group (PDG) average T(a — > KK)/T(ao —> 
7777) = 0.183 ± 0.024 0] to obtain 



/3(a (980) -> 7777) = 0.845 ± 0.017 . 



(4-3) 



For a (1450), we use r(a (1450) -> 777/)/T(a (1450) -»• 777?) = 0.35 ±0.16 and r(a (1450) 
K^)/r(a (1450) -> 7777) = 0.88 ± 0.23 to obtain 



1/B(ao(1450) -»• 7777) = 1.52 ± 0.13 . 



(4.4) 



3 This is in agreement with the value of B(/ o (980) -+ 77+77") = 0.46 ± 0.06 obtained in 0- 
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TABLE IX: Same as Table IIXI except for B decays to a scalar meson and a vector meson. 



Mode 




Theory 


Expt 


Mode 


Theory 


Expt 


B~ -> 


/o(980)A-*- 


o q+0.2+0.7+5.7 
J -°-0.2-0.8-0.9 


-15 ±12 


B U 


->■ / (980)if U 


q ^+0.3+0.1+0.7 
°-°-0.2— 0.1-1.3 


7± 10 


B~ -> 


/o(980)p- 


70 ri+3.2+5.1+5.9 
' J - u -3. 4-7. 0-6.1 




B° 
B° 


-> /o(980)p° 
-> /o(980)w 


£0 9+2.4+0.6+ 22.5 
"°- z -2.3-0.5-138.3 
re n +l. 6+0.5+ 28.5 
Od - u -1.5-0.4-108.4 




B- -> 


ag(980)A"*- 


a K +0.6+0.3+ 2.0 
^• u -0.6-0.3-18.2 




B° 


->■ a+(980)A^*- 


q O+0.6+0.3+ 0.1 
o -°-0.6-0.3-21.6 




B~ -> 


op (980) a: 


n ,1 +0.0+0.0+0.0 
u -^-0.0-0.0-0.3 




B° 


->■ ag (980)1?*° 


n 1 +64.1+72.7+ 2.0 
u - 1 -79.3-73.3-18.2 




B- -> 


ag(980)p- 


in 7 +0.6+l.l + 19.0 
lu - ' -0.7-1.4- 1.4 




B° 


->■ a^(980)p" 


-1 q+0. 1+0.0+8.2 
L - 3 -0. 1-0. 0-0.0 




B~ -> 


_o (980V 


a 9+2.7+1.8+28.1 
^•^-2.1-1.8-12.4 




B° 


-> ap (980)/9+ 


1 9 7+O.7+O.O+ 6.3 
' -0.6-0.0-24.5 




B" -+ 


a„ (980)w 


£ 9+2.9+0.6+47.5 
°- z -2. 1-1.0- 7.1 




B° 
B° 


-> a°(980)p° 
-+ a[](980)w 


7 q+0. 7+0. 2+12. 5 
' ,0 — 0.6-0.1— 11.7 
79 Q+8.9+0.3+23.4 
' z -°-6. 4-0. 3-22.0 




B" -> 


(1450) AT*" 


9 £+0. 6+0.3+ 3.6 
z -°-0.5-0.4-20.2 




B° 


-> a+(1450)A^- 


A0 fi +2.3+4.0+ 0.6 
^•"—2.5-5.1—140.4 




B- -> 


_o (1450)AT 


fl qc+0. 02+0. 02+0. 08 
"■^"-0. 02-0. 02-0. 88 




B° 


-> a8(1450)i_*° 


-1 |-,+0.2+0.0+3.3 
••■•"-0.2-0.1-9.6 




B~ -> 


og(1450)p- 


1 a 1 +1-9+2.7+45.2 
ly - i -2.3-3.8- 5.4 




B° 


-> a+(1450)p- 


1 a +0.1+0.0+19.0 
i '*-0.1-0.0- 0.4 




B" -> 


a„ (1450)p° 


1 e+1.5+1.3+37.2 
J -°-l. 2— 1.3-18.4 




B° 


-> 00 (1450)^+ 


1 q q + 0. 9 + 0.0+ 7.1 

iJ -°-0.8-0.0-22.9 




B" -> 


a (1450)w 


1 fi +2.3+1.5+54.7 
liU -l. 7-1. 7-18.1 




n 

B° 
B° 


->• a[](1450)p 
ag(1450)a; 


£ ^+0.7+0.3+12.2 
"• -0.6-0.3- 5.3 

— Ol. »_2.9-1.0— 21.3 




B" -> 


K*- (1430)^ 


R/l+0. 02+0. 01+0. 32 
"•"^-0.03-0.01-1.49 


4± 15 


B° 


->■ A ? o°(1430)^ 


n a q+0. 04+0. 01+ 3.61 
u -^°-0. 04-0. 00-20. 2 


20 ± 15 


B" -> 


^°(1430)p- 


n 09+0.34+0.33+0.71 

u - oz -0. 29-0. 31 + 1. 09 




B° 


->_f *-(1430)p+ 


1 1+0.0+0.2+0.2 
A,A -0.0-0.2-1.2 




B" -> 


AT*-(1430) j o° 


1 fi +0. 6+0. 1+2.8 
J -°-0.6-0.1-6.2 




B° 


-> i?;°(1430)p° 


n r/i +0.45+0.02+3.76 
u - d ^-0. 46-0. 02-1. 80 




B" -> 


A'*"(1430)cj 


n cr+0.35+0.08+0.32 
v -°°-0. 34-0. 08-2. 49 


-10 ±9 


B° 


-> Zo (1430)a; 


n nq+0. 37+0. 01+0. 29 
U.UO_ . 35-0. 01-3. 00 


-7±9 



A. Decays involving a _vq(1430) state 

Among the hadronic 2-body B decays with a l_g(1430) in the final state, only B — > KqIjO 

BK* BK* 

and B —> K^cj) are sensitive to the B —> iTo(1430) transition form factors F and F 1 °, 
respectively. It turns out that the measurement of B — > K^rfi^ favors a smaller Fq K ° (0), while the 

BK* 

data of B — > Kq4> prefer a smaller F 1 (0) inferred from the neutral mode and a slightly larger 

BK* 

_2 (0) from the charged mode. Since i*i(0) = Fq(0), this means that it is preferable to fit the 
data of B — > K^rf^ and B — > K^cf) by a smaller B — > ^(1430) transition form factor. In this work 



we shall use Fq K ° (0) ~ 0.26 obtained in the covariant light-front quark model |32|. 



1. B -> _fJ(1430)7/C) decays 



Following [8|, we write 

A B ^K*ri' = x( BK °^C l + X { - BK ^C 2 + xW*' K o)Cz, (4.5) 

where C\, C 2 and C3 terms correspond to Figs. 1(a), 1(b) and 1(c), respectively, and X^ BS,P \ 
X(bp,s) are f ac torizable terms defined in Eq. ()A6p . The expressions of Cj's in terms of the 
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s 

~<- 

s 

->— 



(b) 



Us 



K\ 



<1 



FIG. 1: Three different penguin contributions to B 
penguin operators 03,5,7,9. 



K* (U30)rj^\ Fig. 1(a) is induced by the 



parameters can be found in [8]. Because of the small vector decay constant of ii"o(1430), 



x (B v p,K*) is SU pp resse d relative to X^ BK S^P) an d X^ BK o^ ] ). 



However, the C3 term gains a 
= 12.3 at n = 4.2 GeV 



large enhancement from a%{r] q KQ) due to the fact that the chiral factor 
is larger than r/F = 1.5 by one order of magnitude owing to the large mass of -?Tq(1430). It follows 
that a%{r} q KQ) is much greater than a^K^ris) and a^(KQT] q ). As a result, the amplitude of Fig. 
1(c) is comparable to that of Fig. 1(a). 

Because of the large magnitude of a3(K^r] q ^) and the large cancellation between a^K^ris) and 
04(^0*77,5) in C2, B — > KqT]^ decays are dominated by the contributions from Figs. 1(a) and 1(c) 
[8j. Therefore, the penguin diagrams Figs. 1(a) and 1(c) contribute constructively to both K^rj and 
K$rf with comparable magnitudes. Since x( BK o>Vi) /x( BK o ><) = l( B l^S)/lW^) = C ot <p « 
1.23, it is clear that Ab^-k't)/ Ab^k'ti' ~ cot^ and hence B — > K^r] should have a rate larger than 
B —7- KqT]' in scenario 2 as the mixing angle <j) is less than 45°. 

As mentioned in the passing, we use the form factor Fq K ° (0) ~ 0.26 in this work. If a large 
form factor, say, F^ K ° (0) ~ 0.45 H, is employed, we will have B(B -> (1430)r/) ~ 21 x 10~ 6 > 
B{B -)• ^(1430)7/) ~ 3.2 x 10" 6 . This indicates that a small form factor for B 1^(1430) 
transition is preferable in this case. 

A recent pQCD calculation |]J] shows B(B Kqt/) « 75 x 10~ 6 » /3(S ->• A'q??), m sharp 
contrast to the experimental measurements £>(-£> — > KqT]') < B(B — > K^rj) and B(B — > K^r]') ~ 
6 x 10~ 6 (see Table EJ). 



2. B -+ ET *(1430)7r decays 

It is clear from Table IVl that while the predicted branching fraction of B° — > i^o°(1430)7r° is 
consistent with the data, the calculated rates of K (1430)7r~ and -K"o _ (1430)7r + are too small 
compared to experiment. Under the isospin limit, it is naively expected that 

_ B(B° -» Ko°(1430)7r°) _ 1 _ B(B~ -» if*-(1430)7r°) _ 1 

B(B° -> K*-(1430)vr+) ~ 2' 2 ~ ^ 7^ (i430)7r-) ~ 2' 
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T(B ») S (B-^<(1430) T -) =1 
t(B-)B(B° if *-(1430)7r+) 
However, the first two relations are not borne out by experiment. In principle, the data of 
i^g _ (1430)-7r + and K (1430)-7r _ can be accommodated by taking into account the power cor- 
rections due to the non-vanishing pa and pn from weak annihilation and hard spectator inter- 
actions, respectively. However, this will affect the agreement between theory and experiment for 
B -> ^(1430)t? (,) and B i^(1430)p. Indeed, a global fit of p A and 4> A to the B -> SP data 
shown in Table IVl yields pa = 0.15 and 4>a = 82° with x 2 = 8.3. The best fitted results are 
B(B- K* °(U30)ir-) ~ 13.4 x 1(T 6 and B(B° -> K^(U30)tt + ) ~ 14.1 x 1(T 6 which are very 
close to the QCDF predictions with p A = and 4>a = 0. 

The measured branching fractions of B° — > Kq~(1430)7t + and B~ — > i^g°(1430)7r _ are of 
order 30 x 10~ 6 and 50 x 10~ 6 by BaBar and Belle, respectively (see Table |T]), though they are 
consistent with each other within errors. This is probably ascribed to the fact that the definitions 
of the i^g(1430) and nonresonant contributions by BaBar and Belle are different. While Belle 
employed the Brei-Wigner parametrization to describe the -?Tq(1430) resonance, BaBar used the 
LASS parametrization to describe the Kir S-wave and the nonresonant component by a single 
amplitude suggested by the LASS collaboration [44| to describe the scalar amplitude in elastic 
Kir scattering. Since the LASS parametrization is valid up to the Kir invariant mass of order 
1.8 GeV, BaBar introduced a phase-space nonresonant component to describe an excess of signal 
events at higher Kir invariant mass. Hence, the BaBar definition for the 1Tq(1430) includes an 
effective range term to account for the low Kir S-wme while for the Belle parametrization, this 
component is absorbed into the nonresonant piece. In order to compare the BaBar results with 
the Belle ones determined from the Breit-Wigner parametrization, it would be more appropriate to 
consider the Breit-Wigner component only of the LASS parametrization. Indeed, the BaBar results 
for B — > i^o(1430)7r quoted in Table IVl are obtained from (Kir)Q ir and (i^7r)o~vr + by subtracting 
the elastic range term from the Kir S-wave 0, 45]. However, the discrepancy between BaBar and 
Belle for the K^ir modes still remains and it is crucial to resolve this important issue. 

Contrary to the decay mode B — > i^Q _ (1430)7r + whose rate is predicted too small compared to 
the data, the calculated B{W -)• K^~(U30)p + ) is in accordance with experiment. Therefore, it is 
a puzzle why the QCDF approach works well for Xq _ (1430)/9 + but not for i^Q _ (1430)7r + , whereas 
it is the other way around for pQCD. 

It appears that the calculated branching fractions of K ir~ and K^~ir + by pQCD 14J are in 
better agreement with the data (see Table W\ . However, the predicted rates of B — > K^p in this 
approach are too small as we shall see shortly below. 



3. B -+K% (1430) (p,uj,(j)) decays 

We see from Table IVT1 that the calculated B — > Kq (1430) (0, p, u) rates are in good agreement 
with experiment, though the central value of B~ —> Kq~ (1430)0 is smaller than the data. It is 
obvious that the data of B — > Kq (1430) (0, p, oj) can be well accommodated without introducing 
penguin annihilation effects characterized by the parameters pa and 4>a to the central values. 
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Note that the predicted branching fractions of B — > .K"q(143O)0 and B — > i^"g(1430)p in this 
work are substantially smaller than those shown in Table IV of [7|. We found that the rate of 
B — > Kq (1430)0 is sensitive to the scale of // and is large at p = rrib/2 which is the scale used in 
0| . While B{B — > Kq(1A30)p) is stable against the choice of p,, we found a sign mistake in the 
computer code of the previous work [7j; the sign in front of $^(77) in Eqs. (3.7) and (3.8) of [7| 
should read =F instead of ±. As a consequence, the calculated B(B — > Kq(1430)p) in [7|] were too 
large. 

As stated before, we use F (0) ~ 0.26 in order to account for the data of B -»• if *(1430)r/O. 
The predicted B -»• K (1430)0 is in good agreement with experiment, while the central value of 
B(B~ — > Kq~ (1430)0) is smaller than experiment. Nevertheless, we notice that the current data 
imply a large disparity between iTo°(143O)0 and Kq~ (1430)0. 



For comparison, we see from Table IVII that in the pQCD approach 15|, [18|], the predicted 
branching fractions are too small for B —> Kq (1430)(p, 00) and too large for B — > Kq (1430)0 
compared to experiment. It is interesting to see that the rates of B — > KqTt are larger than that of 
B — > KqP within the pQCD framework, whereas it is the other way around in QCDF. 



B. Decays involving a ciq meson 

As stated before, we shall work in scenario 1 for the light scalar meson ao(980) and scenario 2 
for the heavy one ao(1450) so that both of them are ground qq bound states. We see from Table IVl 
that in general B{B — > ao(980)K) is only of order 10~ 7 . 4 This may explain why ao(980) has not 
been seen thus far in hadronic B decays, whereas plenty of ao(980) events have been observed in 
D decays. Notice that B(B — > ao(980)K) are predicted to fall into the range of (4 ~ 10) x 10~ 6 in 
the pQCD approach [12|] and all of them are ruled out by experiment except B — > clqK° '. 

Contrary to B — > PP decays where the production of ir + K~ is substantially greater than 7r + 7r~ , 
it is expected in QCDF that ao(980) + K~ has a rate much smaller than ao(980) + 7r~. Consider 
the interference between the QCD penguin amplitude governed by a^(aoP) — r^a^aoP) and the 
penguin annihilation amplitude proportional to (V u bV* q + V c bV* q )/s/a /p&3(ao-P) for P = K, ir with 
q = d,s. Since 04^, 63 and f a + are all negative, Y^p^p = 0-04 and Y^ p ^p = —0.008 — 0.003z, 
the interference is destructive for 5° — > clqK~ and constructive for 5° — > a f 7r~. This explains 
why the former has a rate much smaller than the latter in QCDF. By contrast, pQCD predicts the 



other way around [l2j, [171] . 



Since the decay amplitudes of the tree-dominated decays B° — > Oq1t~ and B° — > aQ~ir + are 
proportional to aif n F Ba ° and aif a -F B7T , respectively, it is thus expected that B(B° — > Oq-k^) S> 
B(B° — >• <2q 7r + ) for both ao = ao(980) and ao(1450) as the decay constant f a - is very small. 
Moreover, the mode clq-k - should have the largest rate among various B — > aovr decays. Therefore, 



the pQCD results for a^(980)vr- and a£(980)ir~ [1J are at odd with the expectation: The rate of 



In 6], the predicted rates are too large for ao(980)if in scenario 1 and too small for «o(1450)A' in scenario 
2. This is ascribed to a typo appearing in the computer code for the chiral factor rff . 
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ag(980)7r is predicted to be larger than ao~(980)7r and a (980)7r + is not very suppressed relative 
to a+ (980)7r-. 

The decay rate of B° — > a^(980)7r =F is not simply the sum of T(B° — > ag~ (980)71"-) and T(B° —¥ 
Oq (980)7r + ) because of the interference of ao~(— y W^) 77 ^ with (— > 7/7r~)7r + . In principle, one 
should study the 3-body decay B — > t]it + it~ and then apply the narrow-width approximation 

B(B° -> r]ir + ir-) = B(B° -> a^ T )B{a± ? ? 7r ± ) (4.7) 

to determine the branching fraction of B° — > a^ir^. While a direct evaluation of the decay am- 
plitude of B° — > 7]ir + -K~ can be done in the factorization approach (see e.g. 46]), we will simplify 
the calculation by assuming that the amplitude A(B° — > a^7r T ) is the sum of A(B° — > a^7r _ ) and 
A(B° — > ag 7r + ). Since the decay constants f a + and f a - are of opposite sign [see Eq. (|2.4|) ]. the inter- 
ference between clq1t~ and ag tt + is destructive and we find B(B -»■ a (980) ± 7r^) = (4.7+f;|) x lO" 6 
in QCDF. Within theoretical uncertainties this is consistent with the current experimental limit 
3.7 x 10~ 6 set by BaBar [47]. Likewise, the calculated B(B° -> a (1450) ± 7r T ) ~ 1.3 x 10" 6 is also 
in accordance with the limit 3.5 x 10~ 6 [471 ]. 

So far the results of B(B — > ao(980)7r) and B(B — > ao(980)K) in QCDF are all consistent with 
the experimental limits. Hence, at this moment, we cannot conclude on the 2-quark or 4-quark 
nature of ao(980). Nevertheless, as stressed in [6J], if the measured rate of a^(980)7r =F is at the 
level of 1 x 10~ 6 or even smaller, this will imply a substantially smaller B — > ao(980) form factor 
than the B — > ir one. In this case, the four-quark explanation of the ao(980) will be preferred to 
account for the B — > ao(980) form factor suppression. Since ao(1450) can be described by the qq 
quark model, the study of a^(1450)7r =F relative to a^(980)7r =F can provide a more strong test on the 
quark content of ao(980). We see that if the branching fractions of both B° — > a^ : (1450)7r =F and 
B° —7- a^ : (980)7r =F are measured at the level of 1 x 10~ 6 , it will be likely to imply a 2-quark nature 
for do (1450) and a four-quark assignment for oo(980). 

Among the decays B — > ao(980)p and B — > ao(1450)/9, it is clear that a[j~(980)/9~ and aj~(1450)p - 
have the largest rates. Since the branching fraction of B° — > a^(980)p T is predicted to be large 
in QCDF, of order 23 x 10~ 6 , a measurement of this mode may give the first observation of 
the ao(980) production in B decays. Notice that the calculated rates of B~ — > ao(1450)/?~ and 
B — > a^(1450)p _ in the pQCD approach 19(] are unreasonably too large. 



C. Scenario 1 for a (1450) and A" *(1430) 

In Table IVIII we show the branching fractions of B — > SP, SV decays with S = 
ao (1450), Kq (1430) in scenario 1 where scalar mesons ao(1450) and Xq(1430) are treated as the first 
excited states of low lying light qq scalars ao(980) and Xq(800). It is evident that this scenario for 
heavy scalar mesons is ruled out by experiment. 5 For example, the predicted B(B° — > a^(1450)7r =F ) 
is too large and the branching fractions of B — > Kq (1430) (p, u>) are too small compared to the data. 



5 An exception is the decay B — > 1^(1430)77': the calculated rates in scenario 1 are consistent with the 
data. 
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In Tables [V] and IVII we have found a better agreement between theory and experiment for scalar 
mesons above 1 GeV in scenario 2 in which ao(1450) and Kq(1450) are lowest lying qq states. 
This also implies that the light scalars Xq(800), ao(980) and /o(980) are preferred to be four-quark 
bound states. 



D. Decays involving a /o(980) meson 

The penguin- dominated B -)• / (980).K"W decays receive three distinct types of factorizable 
contributions: one from the K^*' emission, one from the /o emission with the ss content, and the 
other from the /o emission with the nn component; see Eq. dH]). Therefore, B{B -)■ / (980)ifM) 
depends on the mixing angle 9 of strange and nonstrange components of the /o(980). In Fig. [2] 
branching fractions of B~ -> / (980) If" and B~ -> / (980)K* _ are plotted versus the /o(980) - 
/o(500) mixing angle defined in Eq. (|2,11|) . We see that both foK~ and foK*~ rates can be 
accommodated with 6 in the vicinity of 20° without introducing 1/rrib power corrections from 
penguin annihilation. This range of the mixing angle is consistent with the limit \0\ < 30° set 
recently by LHCb. 

For definiteness, we show in Tables IVllIXI the branching fractions and CP asymmetries of B — > 
/ (980)irW for 6 = 17°. With this mixing angle, the calculated branching fraction of B — > 
/o(980)p° by QCDF or pQCD is much smaller than the Belle measurement B(B -> / (980)p°) = 
(1.7 ± 0.6) x 10~ 6 as shown in Fig. [3l Since the BaBar measurement yielded only an upper 
bound B(B -»• / (980)p°) < 0.8 x 10" 6 [43], the experimental issue remains to be resolved. 

In order to make quantitative calculations for B — > fo(980)K^*\ we have assumed the conven- 
tional 2-quark description of the light scalar mesons. However, the fact that their rates can be ac- 
commodated in the 2-quark picture for /q(980) does not mean that the measurements of B — > foK^*' 
can be used to distinguish between the 2-quark and 4-quark assignment for /o(980). As discussed 
in 0, 149I ]. the number of the quark diagrams for the penguin contributions to B — > fo(980)K^ 
in the four-quark scheme for /o(980) is two times as many as that in the usual 2-quark picture. 
Therefore, there is no reason that the B — > fo(980)K^ rate will be suppressed if /o is a four-quark 
state. However, in practice, it is difficult to give quantitative predictions based on this scenario 
as the nonfactorizable diagrams are usually not amenable. Moreover, even for the factorizable 
contributions, the calculation of the /o(980) decay constant and its form factors in the four-quark 
scenario is beyond the conventional quark model. 



E. CP violation 



Thus far CP violation has not been observed in any B decays involving a scalar meson. The 
predictions based on QCD factorization are summarized in Tables IVIIII and IIX1 Mixing-induced 



CP asymmetry in B — > fo(980)Ks decays has been studied in [6j, [50j, [51|. 
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FIG. 2: Branching fractions of (a) B~ -> / (980)K" and (b) B -> /o(980)K*" versus the 
mixing angle 6 of strange and nonstrange components of /o(980), where the middle bold solid 
curve inside the allowed region corresponds to the central value. Theoretical errors due to the 
penguin annihilation are taken into account. The horizontal band within the dashed lines shows 
the experimentally allowed region with one sigma error. 




FIG. 3: Same as Fig. [2] except for the decay B° -»■ / (980)p°. 



V. CONCLUSION 

We have studied the hadronic charmless B decays to scalar mesons within the framework of 
QCD factorization. The main results are: 

• We have considered two possible scenarios for the scalar mesons above 1 GeV, depending 
on whether the light scalars Kq(800), ao(980) and /o(980) are treated as the lowest lying 
qq states or four-quark particles. We found that the experimental data favor the scenario in 
which the scalar mesons ao(1450) and -K"q(1430) are the lowest lying qq bound states. This 
in turn implies a preferred four-quark nature for light scalars below 1 GeV. 
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The data of B ^(1430)/) and 5 -> £T *(1430)(p, uj,<p) can be accommodated within 
the framework of QCD factorization without introducing power corrections from penguin 
annihilation, while the predicted B~ — > i^Q°(1430)7r _ and B° — > A'q _ (1430)7t + are too small 
compared to experiment. In view of the fact that the calculated K^p rates are in good 
agreement with experiment, it is very important to have more accurate measurements of 
B —> KqIt decays to pin down the discrepancy between theory and experiment for K^n 
modes. 

If i^g(1430) is made of the lowest-lying qq, we found that Figs. 1(a) and 1(c) interfere 
constructively and that A(B — > Kq1])/A(B — > K^rj') ~ cot</> with (j) being the rj-rj' mixing 
angle in the rj g and rj s flavor basis. Hence, K^r] has a rate slightly larger than K^r/' owing to 
the fact that cj> is less than 45°. This is in sharp contrast to the B — > Kr( decay which has 
the largest rate in 2-body decays of the B meson. 

To accommodate the data of B — >■ i^Q (1430) (77, 77') and B — > Kq (1430)0 we found that a 
small form factor for B — > 1Tq(1430) transition as obtained in the covariant light-front quark 
model is preferable, though other approaches such as pQCD and QCD sum rules tend to 
yield large form factors for B to S transitions. 

We have corrected the results for ao(980)K and ao(1450)iT modes obtained in the previous 
study Branching fractions should be of order 10 -7 for B — > a>o(980)K in scenario 1 and 
of order 10 -6 for B — > ao(1450)-ftT in scenario 2 rather than 10~ 6 and 10~ 7 , respectively, as 
predicted before. It is expected in QCDF that ao(980) + K~ has a rate much smaller than 
ao(980) + 7r~, whereas it is the other way around in pQCD. Experimental measurements of 
these two modes will help discriminate between these two different approaches. 

Although it is widely perceived that light scalar mesons such as /o(980), ao(980) are pre- 
dominately four-quark states, in practice it is difficult to make quantitative predictions on 
B — > SP based on the four-quark picture for S. Hence, in practice we shall assume the two- 
quark scenario for light scalar mesons in calculations. So far the calculated B — > ao(980)-7r 
and B — > ao(980)K rates in QCDF are all consistent with the experimental limits. Hence, we 
cannot conclude on the 2-quark or 4-quark nature of ao(980). Nevertheless, if the branching 
fraction of B° — > a^(980)-7r =F rate is found to be smaller, say, of order 1 x 10~ 6 , it could imply 
a four-quark assignment for ao(980). Since B(B° — > a^(980)p T ) is predicted to be large in 
QCDF, of order 23 x 10~ 6 , a measurement of this mode may give the first observation of the 
ao(980) production in B decays. 

Assuming 2-quark bound states for /o(980) and /o(500), the observed large rates of /o(980)K 
and /o(980)i^* modes can be explained in QCDF with the /o(980)— /o(500) mixing angle 9 in 
the vicinity of 20°. However, this does not necessarily imply that a 4-quark nature for /o(980) 
is ruled out because of extra diagrams contributing to B — > fo(980)K^*\ Irrespective of the 
mixing angle 6, the predicted rate of B — > /o(980)p° is far below the Belle measurement and 
this needs to be clarified in the future. 
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• Contrary to the odd-parity meson sector, we found that penguin-annihilation induced power 
corrections are not needed to explain the penguin dominated B — > SP and B — > SV decays 
in QCDF except for A^(1430)vr modes. How to understand both i^(1430)vr and A *(1430)p 
simultaneously remains an issue in QCD-inspired approaches. 
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Appendix A: Decay amplitudes of B — > SP, SV 



Within the framework of QCD factorization, the decay amplitudes of B — > SP, SV decays can be 
found in [6] and [7], respectively. However, some factorizable terms were missed in the expressions 
for the decay amplitudes of B — > (/o, af])^, 7r) given before in 0]; that is, we did not take into 
account the contributions from the /o or the neutral ao emission induced from the four-quark 
operators other than Oq and Og (see also [12J). They are corrected here: 6 



V * p=u,c 



3,EW 



-Of 



4,EW 



+ SbSk 



A(B^f K") = ^ e ^ s) \ (^+«s- 



ffz(h&u + b 3 + 6 3 ,ew) fu K + ff°{b 2 5 p u + b 3 + &3,ew)jc/« 



V^2 p =u c 



+ 6c 



-a 



3,EW 



— oc 



4,EW 



X 



V 

4,EW 



+ SbSk 



h 
2 



ff$( b 3 ~ 2 fo 3,Ew) fd K + ff*(b 3 - -h,Ew) K f° 



6 In Eq. ()Aip . when applying bar to 0:4 and o^ew, we do not apply bar to ag and as; that is, for those 
Wilson coefficients associated with r x , the bar applies only to r x itself. For example, 0:4=04 =F r x ae- The 



06,8 terms in Eq. (|A1[) are in agreement with those in the formulas in Appendix A of [6|. 
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A(B~ -»• a° K-) = 



G F 



p=u,c { 



+ (a2)Ka X iBK ' ao) + f B fKfa (h6 P u + h + &3,E W ) 



a K ( ) 



p=u,c 



V 

*4,EW 



+ ( a 2 5 p + 2af + -a| EW + a£ - \a 



4,EW 



fBfnffS ( b ^l + h + &3,EW) + (M£ + &3 + &3,Ew) T /j- 



A(S U / 7T° 



p=u,c 



E A ? + «4 " |«3,EW " J«4,EW * ( *# 



+ ( a2<^ + 2a§ + + -a§ jEW - -a 



p 

4,EW 



d 



- - 1- 3- - - 1- 3- 

{biK ~h + 2 & 3,ew + -b 4 ,Ew) f d n + (W£ - &3 + 2 & 3,ew + 2 & 4,ev 



p=u,c 



-(Bn,ao) 



A(B° -> a°vr ) = E/fj - «S + ^< EW + ^ew)^^ '* 



+ (^a 2 <5£ - a p 4 + 2<,ew + 2 a Wj X 

+ f B fJa [(hSP + h + h,Ew) aon ~ (hS? + h + &3,Ew) wao 



f ( a 2 <S£ - a 4 + t^ew + 



-^(B7r,a ) 



7rao 



fBfwfao (M£ + &3 + 26 4 - ~(&3,EW ~ &4,E\v)) 



ao7r 



+ (M£ + &3 + 26 4 - tj(&3,ew - Hew)) 



7rao 



Here ao stands for ao(980) or ao(1450), X p = V pb V* q with q = d,s 



r?(M) 



2m^ 



K* 



2mi. 



m b (n)(m s (n) - m q (n)) ' 



2m 



CI(> 



x m b (fi)(m d (iJ,) - m„ (//))' ' x 



f a °(u) = 2mao f fo (u) = 2mfo 
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and 



c^(MiM 2 ) 



a p z (MiM 2 ) + of (MiM 2 ); for MiM 2 = PS, SV, VS, 
4 (MiM 2 ) - af (MiM 2 ); for MiM 2 = SP, 
<(MiM 2 ) - rf 2 ag(MiM 2 ); for M X M 2 = PS,SP,SV , 



4V 1 z / I „P T\,f Ti/r \ 1 „Mo „V ~MT A/T \. t I\/T )\/T T^C <■ ' 



«3,EW( M 1 M 2) 



og(MiM2)+r^ 3 ag(MiM 2 ); for M l M 2 = VS 
al(M x M 2 ) + a p 7 {M 1 M 2 ); for M X M 2 = PS, 5V, VS, 
al(M 1 M 2 ) - a p 7 (M 1 M 2 ); for M X M 2 = SP, 



, „ „ , f a? (MiM 2 ) - rf^KMiM,); for MxM 2 = PS, SP, SV, 
\< (M 1 M 2 ) + rf 2 af(M 1 M 2 ); for M l M 2 = VS. 



Q P EW (M 1 M 2 ) 

The general expression of the effective Wilson coefficients aj is given in Eq. (|3.ip . while the 
effective Wilson coefficients 04 appearing in Eq. (]Aip are defined as aifig 1 and they can be obtained 
from Eq. (|3.ip by retaining only those terms that are proportional to ^5. Specifically, 

a^M.M,) = ^ ^ [V(M 2 ) + ^^( Ml M 2 )] + Pf (M 2 ). (A4) 

The LCDA of the neutral scalar meson is replaced by $5 which has the similar expression as Eq. 
(|2.8|) except that the first constant term does not contribute and the term f$ us is factored out: 

00 

$ s (x,ii) = 6x(l-x) BMCU 2 {2x-l). (A5) 

m=l 

The annihilation terms bi have the same expressions as bi except that r'^ and u-sBi are replaced by 

^dS^pec^Mseea). 

In Eq. (|A1|) . the quantity X is a factorizable term whose explicit expression is given by 



2 1 



x (bs,v) = (y|jM| 0)( 5|^|^ )=2 / ymi3 p c F 1 B5 (m : 

= (5|J"|0)(V|^|B) = -2f s m B pX V (rn%), (A6) 
X (BS,P) s (P |jM|o)(S|j;|P) = -/p(m|-m|)P ( f s (m 2 p), 

= (5 | jM|o)(P|j£|B) = / s ( m | - m P )P ^(m|), 

where p c is the cm. momentum. 

Note that the /o(980)-/o(500) mixing angle (i.e. sin#) and Clebsch-Gordon coefficient l/y/2 
have been included in the /o(980) form factors F B ?o and decay constants /j£ • Throughout, the 
order of the arguments of the aF i (M\M 2 ) and bi(M\M 2 ) coefficients is dictated by the subscript 
M\M 2 , where M 2 is the emitted meson and M\ shares the same spectator quark with the B meson. 
For the annihilation diagram, M\ is referred to the one containing an antiquark from the weak 
vertex, while M 2 contains a quark from the weak vertex. 
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